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一 种 考虑 时 变温 度 场 影 响 的 非 线 性 闸 振 计算 方法 


沈 恩 楠 ,吴江 鹏 , 曲 林 锋 , 王 赫 直 , 赵 飞 宇 
(中 国航 空 工业 集团 公司 沈阳 飞机 设计 人 研究 所 ,110035 沈阳 ) 


摘 要 :基于 CFD 的 颤 振 计 算 时 域 方 法 通常 在 模 态 空间 内 耦合 求解 Euler/N-S 方程 和 结构 运动 方 
程 ,局 限于 线性 结构 。 为 了 研究 温度 等 非 线 性 因素 对 颤 振 特性 的 影响 ,本 研究 将 Nastran 软件 的 线 
性 和 非 线性 结构 运动 方程 逐步 积分 求解 的 CSD 方法 与 非 定常 CFD 求解 器 时 域 耦 合 , 并 利用 径 向 基 
函数 法 进行 流 国 耦合 界面 数据 传递 ,构建 出 结构 线性 和 非 线性 的 颤 振 计 算 方法 。 采 用 所 建立 的 结 
构 线 性 和 非 线性 颤 振 计算 方法 对 一 平板 机 翼 进 行 分 析 , 验 证 了 方法 的 可 靠 性 。 最 后 建立 了 能 够 模 
拟 温 度 实时 变化 过 程 的 颤 振 计算 方法 ,研究 了 温度 对 材料 的 影响 引起 的 非 线性 颤 振 特性 的 变化 。 
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Nonlinear flutter calculation method considering the 
influence of time-varying temperature field 


SHEN Ennan, WU Jiangpeng, QU Linfeng, WANG Hezhe,ZHAO Feiyu 


(Shenyang Aircraft Design and Research Institute of AVIC,110035 Shenyang, China) 


Abstract: The time domain method of flutter calculation based on CFD usually solves Euler/ N-S equations 
and structural motion equations in modal space ,which is limited to linear structures. In order to study the 
influence of temperature and other nonlinear factors on flutter characteristics ,in this paper, the CSD meth- 
od of solving the linear and nonlinear structural motion equations of Nastran software is coupled with the 
unsteady CFD solver in time domain. The radial basis function method is used to transmit the fluid struc- 
ture coupling interface data. The structural linear and nonlinear flutter calculation methods are construc- 
ted. The linear and nonlinear flutter calculation methods are used to analyze a flat swept wing. The results 
verify the reliability of the methods. Finally ,the flutter calculation method which can simulate the real-time 
temperature change process is established, and the non-linear flutter characteristic change caused by the 
influence of temperature on the material is studied. 


Key words flutter; nonlinearity ;time-varying temperature field; coupling 


闸 振 预测 是 飞行 器 设计 中 必 不 可 少 的 一 项 重要 ( computational fluid dynamics ,CFD ) 和 计算 结构 动力 
环节 。 国 外 在 颤 振 数值 计算 方面 已 经 做 了 大 量 研 学 ( computational structural dynamics ,CSD ) 最 新 成 就 
究 “, 现 阶段 主要 研究 特点 是 采用 计算 流体 力学 ”耦合 求解 非 线 性 尤其 是 气动 非 线 性 颜 振 问题 。 
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国内 近 些 年 来 发 展 的 CFD/CSD 耦合 计算 方法 
中 ,多 数 将 气动 力 模型 或 者 结构 模型 进行 线 化 处 理 。 
刘 南 等 ”基于 频 域 颜 振 分 析 V-g 方法 提出 了 一 种 快 
速 确定 参数 变化 影响 的 高 效 颤 振 预 测 分 析 方 法 ,并 
应 用 于 跨 声速 颤 振 边 界 分 析 。 安 效 民 等 对 二 阶 时 
间 精 度 的 CFDZCSD 耦合 算法 进行 了 研究 ,在 保证 计 
算 精 度 的 基础 上 明显 提高 了 计算 效率 ,并 保持 了 松 
耦合 方法 模块 化 的 特点 。 叶 正 寅 等 "9 利用 双 时 间 推 
进 方法 求解 三 维 非 定常 N-S 方程 ,得 到 与 任意 非 定 
常 运动 对 应 的 非 线性 气动 力 ,采用 模 态 释 加 法 描述 
结构 振动 变形 , 计算 了 风 洞 试验 模型 的 颤 振 特性 。 
陆 志良 等 采用 弹性 变形 技术 生成 块 结构 动态 网 
格 ,CFD 计算 非 定常 气动 力 , 模 态 关 加 法 描述 结构 振 
动 变 形 , 并 提出 变质 量变 刚度 的 方法 计算 飞机 跨 声 
速 题 振 特 性 。 张 苏 华 等 ”采用 CFD/CSD 耦合 方法 
证 算 了 AGARD 445.6 机 翼 的 跨 声 速 闸 振 特性 ,利用 
ANSYS 软件 的 MFX 多 场 耦合 求解 器 求解 非 定常 气 
动力 ,采用 直接 数值 积分 法 获得 结构 响应 。 仲 继 泽 
等 习 基 于 弹性 体 动 网 格 技术 ,发 展 了 一 种 用 于 机 辟 
浇 场 网 格 变形 的 降 阶 算法 。 将 流 场 网 格 所 包围 的 空 
闻 区 域 视 为 虚拟 弹性 休 。 在 保证 计算 精度 的 情况 下 
龙 江 节省 了 计算 时 间 。 侯 良 学 等 "采用 CFD/CSD 
耦合 的 非 线性 方法 对 跨 声 速 颜 振 边 界 进行 预测 , 预 
测 结 果 与 试验 结果 吻合 较 好 。 钮 次 斌 等 "引入 微 
分 环 积 法 ,分 析 高 速 小 展 弦 比 机 缀 的 气动 弹性 问题 ， 
将 济 展 弦 比 机 机 简 化 为 悬臂 板 , 根 据 频率 重合 理论 
对 匣 振 问题 进行 求解 。 夏 健 等 ("采用 流 固 厢 合 方 
法 再 跨 声速 额 振 问 题 进 行 了 数值 模拟 ,在 每 一 个 物 
理 时 间 步 都 采用 聚合 多 重 网 格 方法 ,交替 迭代 求解 
流体 和 结构 控制 方程 。 

综 上 所 述 , 目 前 大 部 分 闸 振 计算 方法 都 是 针对 
线性 结构 ,通过 模 态 玛 加 法 来 描述 结构 振动 变形 ,无 
法 考虑 非 线性 结构 问题 。 本 研究 在 线性 和 非 线性 结 
构 瞬 态 响应 的 逐步 积分 计算 方法 的 基础 上 ,建立 了 
结构 线性 和 非 线性 的 CFDZCSD 时 域 耦合 的 颜 振 计 
算 方 法 ,并 实现 了 能 够 考虑 温度 实时 变化 过 程 的 颤 
振 特 性 分 析 。 


1 结构 瞬 态 响应 逐步 积分 求解 方法 


1.1 线性 结构 瞬 态 响应 
不 计 结构 阻尼 ,在 外 部 激励 的 作用 下 ,结构 运动 
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方程 形式 为 

M u(t) + Ku(t) =P(t) (1) 
其 中 ;M 是 质量 矩阵 ;K 是 刚度 矩阵 ; u (1) 是 加 速度 
问 量 ;u(7) 是 位 移 疝 量 ;P(Y) 是 开 集 向量 。 


对 于 线性 结构 分 析 ,质量 矩阵 M 与 刚度 矩阵 


保持 不 变 。 加 速度 的 差分 离散 表达 式 为 
jn -二 (2) 


其 中 :上 标 n 代表 时 间 层 ;At 是 时 间 步 长 。 

将 式 (2) 代 入 式 (1) 可 得 

2 -2uU" + 

Ar 

在 Nastran 软件 中 ,外 部 激励 己 采 用 相 邻 3 个 时 

间 点 激励 的 平均 值 表示 ,弹性 力也 应 采用 相 邻 3 个 
时 间 点 的 平均 值 。 最 后 整理 为 

Au'' =A,+Au" +Au"™ (4) 


+Ku’=P" (3) 


其 中 


2M K M K 
Ee 

该 方法 类 似 于 Newmark-Beta 方法 。 式 (4) 左 侧 
只 包含 w”, 布 侧 只 包含 w” 和 uw", 通过 迭代 方法 
可 以 求 得 w”*。 

在 每 一 个 物理 时 间 步 ,通过 修改 bdf 文件 的 载 
荷 工 况 卡 片 可 以 用 于 模拟 时 变 载 荷 。 对 于 线性 结构 
的 具体 计算 流程 如 下 。 

1) 利 用 Patran 建立 有 限 元 模型 ,生成 . bdf 文件 。 

2) 将 随时 间 变 化 的 载荷 离散 成 看 干 物理 时 间 
步 , 调 用 Nastran 软件 计算 在 第 一 个 物理 时 间 步 载荷 
作用 下 的 结构 响应 。 

3) 读 取 结 构 响 应 ,生成 新 的 bdf 文件 ,更 改 载荷 
和 边界 条 件 。 

4) 调 用 Nastran 软件 计算 结构 在 新 的 物理 时 间 
步 内 的 结构 响应 。 

5 ) 重 复 第 3) ~4) 步 ,直至 完成 全 部 物理 时 间 步 
的 迭代 步 数 。 


1.2 非 线性 结构 瞬 态 响应 


对 于 非 线 性 结构 ,在 每 一 个 物理 时 间 步 内 ,质量 
和 矩阵 和 刚度 矩阵 会 发 生变 化 ,需要 将 每 个 物理 时 间 
步 内 的 载 匣 问 量 分 解 为 厂 干 段 载 集 增 量 ,逐步 加 载 
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至 结构 的 同时 对 结构 的 刚度 矩 阵 K 进行 更 新 。 具 体 
步 又 如 下 : 

1 ) 利用 Patran 建立 有 限 元 模型 ,生成 . bdf 文件 。 

2) 将 随时 间 变 化 的 载荷 离散 成 知 干 物理 时 间 
步 , 调 用 Nastran 软件 计算 在 第 一 个 物理 时 间 步 载荷 
作用 下 的 结构 响应 。 此 时 载荷 施加 过 程 按 照 如 下 步 
又 进行 。 

(a) 先 施加 总 载 答 的 一 部 分 载 丛 。 

(b) 提取 刚度 矩阵 天。 

(c) 求 解 刚 度 和 矩阵 和 残 值 向 量 ,判断 是 否 收敛 。 

(d) 如 果 收 敛 则 得 到 结构 的 位 移 和 应 力 ; 如果 
不 收敛 , 则 重复 (b) ~ (ce) 步 ,直至 收敛 。 

(e) 继 续 增 加 载荷 增 量 ,重复 (a) ~ (d) 步 ,直至 
100% 载荷 加 载 。 

=3 ) 读 取 结 构 响 应 ,生成 新 的 . bdf 文件 ,更 改 载荷 

利 访 界 条 件 。 

4) 调 用 Nastran 软件 ,重复 第 2) 步 里 的 步骤 , 计 
算 结构 在 新 的 物理 时 间 步 内 的 结构 响应 。 

5 ) 重 复 第 3) ~4) 步 ,直至 完成 全 部 物理 时 间 步 
的 大 代步 数 。 


2CJ 流 - 固 耦 合 边界 信息 传递 


CN 大 多 数 倩 况 下 , 流 场 与 结构 计算 的 空间 离散 方 
法 有 很 大 区 别 。 在 机 履 或 者 舵 面 的 前 缘 、 后 缘 和 冀 
梢 等 外 形 曲 率 较 大 的 地 方 , 流 场 网 格 分 布 比 结构 单 
元 蝎 加 密集 ,需要 在 边界 处 通过 插值 方法 传递 气动 
载 入 和 结构 表面 变形 。 采 用 径 癌 基 消 数 法 (RBFs) 
实 魏 结 构 表面 变形 的 传递 过 程 ,将 结构 单元 点 的 位 
移 畏 值 到 流 场 物 面 网 格 点 。RBFs 构造 的 择 值 函数 
是 一 系列 基 消 数 的 全 加 , 即 


f(x) = py aip (x -x,l) + yx) (6) 


非 定常 流 场 
蕉 进 


结构 变形 


结构 变形 def， 
速度 w， 加 速度 
A 


瞬 态 响应 分 析 


结构 初始 外 形 
det, 
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其 中 :f(x) 为 构造 的 插值 函数 ;n 为 控制 点 个 数 ;a 为 
通过 控制 点 函数 及 附件 条 件 获 得 的 系数 。x, = 
[x ,yw ;zwj] 为 已 知 位 移 的 控制 点 的 坐标 。 0 为 基于 
欧式 范 数 的 基 函 数 ,本 研究 采用 ||x| 作 为 基 瑞 数 。y 
(xX) =66+bix+ by +b3z 为 线性 多 项 式 。 其 中 a,， 
bo ,0 ,bs,b; 为 通过 控制 点 变形 及 附加 条 件 获 得 的 系 
数 。 在 计算 过 程 中 ,虽然 结构 表面 点 的 位 移 在 时 刻 
发 生变 化 ,但 是 与 气动 网 格物 面 点 之 间 的 插值 关系 
是 保持 不 变 的 ,这 样 在 计算 开始 前 就 可 以 确定 并 储 
存 插值 函数 。 在 后 续 的 计算 中 ,结构 发 生 任 意 变 形 
都 可 以 通过 储存 的 插值 函数 快速 计算 出 流 场 物 面 网 
格 的 变形 。 由 于 流 场 网 格 比 结构 单元 分 布 更 密集 ， 
因此 采用 最 近 距 离 法 搜索 流 场 物 面 边界 气动 载 集 到 
结构 单元 表面 之 间 的 对 应 关系 。 


3 ”CFD/CSD 耦合 的 颤 振 计算 


对 于 飞行 器 颤 振 计算 ,外 部 激励 为 气动 表面 压 
强 。 结 构 在 气动 载荷 的 作用 下 产生 变形 , 反 过 来 对 
表面 气动 载荷 分 布 产 生 影 响 。 因 此 ,该 问题 中 的 激 
励 己 是 随时 间 变 化 的 ,时 域 颜 振 计 算 在 时 域内 耦合 
求解 流 场 控制 方程 和 结构 运动 方程 ,根据 结构 响应 
时 间 历 程 判 断 是 否 颤 振 。 本 研究 建立 了 如 下 基于 
CFDZCSD 时 域 耦合 的 颤 振 计 算 方 法 。 

在 气动 力 的 作用 下 ,结构 单元 点 在 每 个 物理 时 
间 ( 耦 合 时 间 At. ) 内 都 会 产生 位 移 .速度 和 加 速度 ， 
并 对 下 一 个 时 间 步 的 结构 运动 方程 求解 产生 影响 ， 
如 图 1 所 示 , 在 每 一 个 物理 时 间 步 (i, ) 求 解 结构 运 
动 方程 时 ,都 要 将 前 一 步 的 计算 结果 作为 初始 条 件 。 
在 求解 线性 或 非 线 性 结构 运动 方程 的 基础 上 ,耦合 
CFD 方法 求解 流 场 控制 方程 ,时 域内 迭代 计算 得 到 
结构 瞬 态 响应 ,进而 分 析 获 得 颤 振 特性 。 


结构 变形 def,， 
速度 ww， 加 速度 


a. 


结构 变形 def,， 
速度 w”， 加 速度 
A, 


图 1 CFD/ACSD 时 域 耦合 颤 振 计 算 流 程 
Fig. 1 Process of CFD/CSD time domain coupled flutter calculation 


投稿 网 站 :http://cjam. xjtu. edu. cn” 微 信 公众 号 :应 用 力学 学 报 


第 6 期 


4 算 例 与 分 析 


首先 以 平板 机 驾 颤 振 分 析 为 例 , 验 证 本 人 研究 所 
建立 的 CFDZCSD 时 域 看 合 的 闸 振 计算 方法 的 正确 
性 ,然后 通过 改变 结构 温度 场 实现 时 域 温度 变化 的 
非 线性 颤 振 计 算 。 


4.1 平板 机 桶 模型 颤 振 分 析 


平板 机 辟 ” 弹性 模 量 为 7. 1 x 10"Pa, 泊 松 比 
0.3 ,密度 为 2700 kg/m  。 翼 根 固 支 ,采用 实体 单元 
建 模 , 展 向 划分 25 个 单元 , 苞 向 划分 16 个 单元 。 前 
后 缘 、 翼 梢 以 及 翼 根 处 网 格 较 密 ,增长 比 为 1.2。 厚 
度 方向 划分 2 个 单元 。 前 后 缘 为 三 棱柱 Wedge6 单 

中 间 部 分 采用 六 面体 Hex8 单元 。 有 限 元 模型 如 
图 多 所 示 。 
Nelli 


有 

©O 

三 

CN 

CN 

CN 

©O 

pi 

> 图 2 平板 机 可 有 限 元 模型 

< Fig.2 Finite element model of plate swept wing 


(6 首先 验证 本 研究 所 建 实体 模型 与 文献 模型 的 结 
构思 动 特性 是 否 一 致 。 利 用 Nastran 软件 分 析 模 态 ， 
前 男 阶 模 态 的 振 型 与 频率 ,如 图 3 所 示 。 通 过 对 比 
( 狗 1) 可 以 看 出 本 研究 所 建立 模型 各 阶 模 态 的 频 束 
与 文献 [ 13] 基 本 吻合 。 


(c) Fre3=278.0 Hz 


(d) Fre4=694.4 Hz 


图 3 平板 机 波 前 四 阶 模 态 


Fig.3 First four modes of the wing 
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表 1 ， 所 建立 实体 模型 与 文献 [ 13 ] 频率 对 比 


Tab.1 Comparison between the presented 


model and Ref. [13 ] 


阶 数 本 研究 模型 /Hz 


文献 模型 /Hz 
i 40.8 39.1 
= 243.2 237.1 
= 278.0 270.5 
四 694.4 704.7 


基于 偶 极 子 格 网 法 建立 非 定 常 气动 力 模型 , 气 
动 网 格 展 向 和 弱 向 分 别 均匀 划分 6 个 和 4 个 单元 ， 
共 24 个 面 元 。 采 用 p-k 法 计算 模型 的 颤 振 速 度 , 空 
气 密度 值 为 1.226 kg/m ,马赫 数 为 0.45。 图 4 所 示 
为 V-g 曲线 ,机 器 的 颤 振 速 度 为 153. 683 m/s, 此 时 
得 到 的 颜 振 速度 是 在 该 来 流 条 件 下 的 匹配 闸 振 速 
度 ,满足 太 =Wa。。 与 文献 [13 ] 计算 结果 ( 颜 振 速 
度 为 149.1 m/s) 吻合 较 好 ,该 有 限 元 模型 能 够 用 于 


颤 振 计 算 方 法 验证 。 
0.4 
02 Ve=153.683 m/s 


阻尼 


140 145 150 155 160 
速度 /(m，s ') 


图 4 ”p-k 法 计算 结果 绘制 V-g 曲线 
Fig.4 V-g diagram from the p-k method 
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图 5 平板 机 可 流 场 网 格 
Fig.3 Conputation mesh of the wing 


采用 结构 线性 的 颤 振 计 算 方 法 分 别 计算 平板 机 


投稿 网 站 :http://ejam. xjtu. edu. cn” 微 信 公 众 号 :应 用 力学 学 报 


1018 应 用 力学 学 报 


辟 在 来 流 马赫 数 为 0.44.0.445 和 0.45 的 情况 。 流 
场 网 格 总 数 153 万 ,对 称 面 及 机 辟 表 面 网 格 如 网 5 所 
示 。 监 测 翼 梢 前 后 缘 顶 点 的 2 方向 位 移 。 当 来 流 马 
赫 数 为 0.44 时 ,监测 点 位 移 曲线 如 图 6(a) 所 示 ; 随 着 
时 间 推 进 ,机 器 振动 幅度 逐渐 减 小 ,监测 点 位 移 曲 线 
收敛 到 静 平 衡 位 置 ; 当 来 流 马赫 数 为 0.445 时 ,监测 


0.10 


位 移 /mm 
位 移 /mm 


(a) Ma=0.44 


思 图 7 为 结构 非 线 性 的 颤 振 计算 方法 得 到 的 来 流 
马 医 数 分 别 为 0.44 .0.45$ .0.453 和 0. 46 的 乙 梢 后 
组 禄 移 随时 间 响 应 昌 线 ,其 中 图 7(a) 与 图 7(b) 包 
会 蝇 线 性 方法 的 计算 结果 。 当 来 流 马赫 数 为 0. 44 
是 淹 线 性 方法 监测 点 位 移 时 间 响 应 曲线 收敛 速度 
明 最 快 于 线性 方法 ; 当 来 流 马 赫 数 为 0.45 时 ,线性 
方法 监测 点 曲线 发 散 ,而 非 线 性 方法 曲线 仍然 处 于 
收 敏 状态 ; 当 来 流 马赫 数 为 0. 455 时 , 非 线性 方法 监 
测 上 点 曲线 已 有 发 散 趋 势 , 颤 振 临 界 马赫 数 在 0.45 ~ 
0. 枢 5 之 间 ; 当 来 流 马 赫 数 为 0. 46 时 , 非 线性 方法 
豚 测 点 曲线 已 经 发 散 , 结 构 发 生 闸 振 。 可 以 认为 非 
线 竹 方法 计算 得 到 颜 振 临 界 马 赫 数 为 0. 453 左右 。 
两 种 方法 计算 结果 基本 一 致 ,其 中 的 差异 主要 来 自 
于 :在 非 线 性 方法 计算 过 程 中 ,载荷 被 分 段 加 载 , 同 
时 刚度 矩阵 在 时 刻 更 新 ,这 种 分 段 积 分 的 形式 与 线 
性 方法 计算 得 到 的 结果 存在 一 定 差 异 。 由 于 平板 机 
翼 颤 振 在 微 幅 振动 条 件 下 不 存在 明显 的 非 线 性 环 
节 , 这 种 差异 并 不 明显 ,所 得 结论 与 物理 客观 相符 。 


4.2 考虑 时 变温 度 影响 的 非 线 性 颤 振 分 析 


当 机 避 结构 温度 升 高 时 ,材料 的 非 线性 就 会 体 
现 出 来 。 主 要 表现 为 材料 的 弹性 模 量 以 及 线 膨胀 系 
数 会 随 温度 的 变化 发 生 改变 ,结构 受热 不 均匀 导致 
各 单元 膨胀 程度 不 同 , 引 发 结构 内 部 的 预 拉 或 预 压 
应 力 ,使 结构 发 生 软 化 或 硬化 。 

结构 线性 的 颤 振 计算 方法 得 到 的 平板 机 疲 海 平 


辟 梢 后 缘 


(b) Ma=0.445 
图 6 结构 线性 的 颤 振 计 算 方 法 得 到 的 惨 梢 前 后 缘 Z 方 向 位 移 曲线 
Fig.6 Z-direction displacement of detection point by CFD/linear CSD coupling analysis 
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点 位 移 曲 线 如 图 6(b) 所 示 , 收 敛 速度 相 比 于 马赫 数 
为 0.44 时 的 情况 减缓 许多 ,接近 赢 振 临界 状态 ; 当 来 
流 马赫 数 为 0.45 时 ,监测 点 位 移 时 间 曲 线 如 图 6(¢) 
所 示 ,结构 振 动 幅 度 逐 渐 增 加 ,结构 发 生 颤 振 。 可 以 
认为 采用 结构 线性 的 颤 振 计算 方法 得 到 的 颤 振 临 界 
马赫 数 约 为 0.445 , 即 闸 振 速度 为 151.3 m/s。 


出 梢 前 缘 ， 站 | 


Wl 


ed 观 梢 后 乡 


位 移 /mm 
© 
[se 


0.4 0.0 0.05 0.10 0.13 
时 间 /s 
(¢) Ma-0 .45 


面额 振 I 临 界 速度 为 151.3 m/s(0.445 Ma) 。 当 来 流 
马赫 数 为 0.4 时 ,结构 不 会 发 生 颤 振 ,此 时 的 来 流速 
压 为 颜 振 临 界 速 压 的 0.808 倍 。 温 度 变 化 可 以 改变 
材料 的 弹性 模 量 ,进而 改变 结构 的 动力 学 特性 , 额 振 
临界 速 压 也 会 发 生变 化 。 所 采用 的 材料 属性 如 表 2 
所 示 。 分 别 计算 结构 温度 为 300 下 .320 K 和 从 300 K 
开始 以 200 K/s 的 速度 均匀 增加 的 情况 。 如 图 8 所 
示 , 当 温度 为 300 K 时 ,监测 点 位 移 曲 线 处 于 收敛 状 
态 ; 当 温度 为 320 K 时 ,位 移 曲线 处 于 发 散 状 态 , 结 
构 发 生 颤 振 ; 变 温 条 件 下 , 当 温 度 较 低 时 ,位 移 曲 线 
处 于 收敛 状态 ; 当 温 度 升 高 至 312 K(0.06s) ~314K 
(0.07 s) 之 间 ,位 移 曲 线 近似 为 等 幅 振荡 ; 当 温 度 升 
高 至 约 为 316 K(0.08 s) 时 ,位 移 曲 线 振幅 明显 增 
加 ,结构 发 生 颤 振 。 由 此 可 以 判断 ,结构 在 312 K ~ 
314 民 之 间 处 于 临界 赢 振 状 态 。 为 了 能 够 更 加 显 性 
地 说 明 温度 对 结构 刚度 的 影响 ,对 320 kK 时 的 结构 
动力 学 特性 进行 分 析 , 前 四 阶 模 态 如 图 9 所 示 。 通 
过 与 结构 温度 300 K 的 情况 对 比 可 以 看 出 ,在 温度 
场 作用 下 ,结构 的 各 阶 模 态 的 振 型 基本 不 变 ,而 频率 
均 降 低 10 儿 左右 ,如 表 3 所 示 。 
表 2 ”材料 属性 随 温度 变化 


Tab.2 Material property varies with temperature 


温度 /KK 300 350 


弹性 模 量 /Pa 7.1x1010 4.0 x1010 
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表 3 频率 变化 对 比 表 


Tab.3 Frequency comparison in different temperature 


阶 数 ”300 K 时 结构 频率 /Hz 320 K 时 结构 频率 /Hz ”变化 /% 
1 40.8 36.9 9.56 
2 243.2 219.8 9.62 
3 278.0 251.3 9.60 
4 694.4 627.6 9.62 
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结构 非 线性 的 颤 振 计算 方法 得 到 的 辟 梢 后 缘 2 方向 位 移 曲 线 
Fig.7 Z-direction displacement of detection point by CFD/nolinear CSD coupling analysis 
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(c) 温度 以 200 K/s 速率 从 300 K 增长 至 320 
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(b) 温度 为 320 KK 

图 8 不 同 温度 状态 下 监测 点 Z 方 向 响应 曲线 


Fig.8 Z-direction displacement of detection point in different temperature conditions 


图 10 ~12 给 出 了 不 同 结构 温度 场 蘑 时 刻 的 流 
线 图 及 温度 和 压强 分 布 图 。 通 过 对 比 图 10 和 图 11 
可 以 看 出 ,由 于 采用 了 等 温 壁 物 面 边 界 条 件 ,结构 温 
度 的 升 高 对 流 场 温度 场 有 一 定 影响 ,但 是 结构 温度 
的 变化 对 流 场 的 温度 场 及 压力 场 分 布 影响 并 不 明 
显 ,温度 场 分 布 的 变化 主要 是 由 结构 较 大 幅度 振动 
引起 的 。 图 12 给 出 了 结构 温度 为 320 K 时 , 某 时 刻 
截面 流 线 图 及 截面 处 的 压力 分 布 , 当 在 发 生 闸 振 时 ， 
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Im 力学 学 报 Cne VD we 


流 场 出 现 较 大 的 分 离 流动 。 对 比 图 13 和 图 14 可 以 
看 出 ,发 生 闸 振 的 机 如 表 面 出 现 大 面积 的 负 压 区 , 属 
于 典型 的 单 阶 模 态 失速 颤 振 。 


(a) Frel=36.87 Hz (b) Fre2=219.81 Hz 


图 12 结构 温度 为 320 kK 时 , 某 时 刻 
流 场 某 截面 流 线 图 及 压力 分 布 


(c) Fre3=251.26 Hz (d) Fre4=627.62 Hz Fig.12 The streamline and pressure distribution 
= 图 9 温度 为 320 K 时 平板 机 器 前 四 阶 模 态 when the structure temperature is 320 K 
Figs9 First four modes of the wing when temperature is 320 K 
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图 13 结构 温度 为 300 kK 时 , 某 时 刻 
图 10 结构 温度 为 300 时 , 某 时 刻 流 场 机 可 表 面 及 附近 压力 分 布 
2 某 截 面 流 线 图 及 温度 分 布 


CFig.10 The streamline and temperature distribution 


naXiv 


Fig. 13 The pressure distribution 


when the structure temperature is 300 K 


C 


when the structure temperature is 300 K 


图 11 结构 温度 为 320 K 时 , 某 时 刻 流 场 
某 截面 流 线 图 及 温度 分 布 机 小 表面 及 附近 压力 分 布 
Fig.11 The streamline and temperature distribution Fig.14 The pressure distribution 
when the structure temperature is 320 K when the structure temperature is 320 K 
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5 结 论 


本 研究 针对 颤 振 问题 建立 了 CFD/CSD 耦合 的 
结构 线性 和 非 线性 颤 振 计 算 方 法 。 以 平板 机 六 模 型 
验证 了 方法 的 可 靠 性 。 进 一 步 采 用 所 建立 的 结构 非 
线性 赢 振 计算 方法 对 时 变温 度 影响 下 的 颤 振 特性 进 
行 分析 , 研 究 了 材料 非 线 性 对 额 振 特性 的 影响 。 得 
出 以 下 结论 。 

1) 采 用 结构 线性 和 非 线 性 的 颤 振 计 算 方 法 得 到 
的 平板 机 可 颜 振 速度 差别 并 不 明显 ,这 与 平板 机 性 
结构 不 存在 非 线性 环节 的 客观 物理 相符 。 本 研究 建 
立 的 颤 振 计算 方法 可 以 用 于 颜 振 预测 。 

2) 温 度 通过 影响 材料 弹性 模 量 对 机 可 刚 度 产生 
影响 。 变 温 条 件 下 , 当 温 度 较 低 时 ,位移 曲线 处 于 收 
敛 状 态 ; 随 着 温度 升 高 ,位 移 曲 线 振幅 明显 增加 , 结 
梅 发生 闸 振 。 结 构 温度 的 变化 对 流 场 的 温度 场 及 压 
为 声 分 布 影响 并 不 明显 ,温度 场 分 布 的 变化 主要 是 
责 厂 构 较 大 幅度 振动 引起 的 。 
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